A few Smarandache Integer Sequences 
Henry Ibstedt 


Abstract 
This paper deals with the analysis of a few Smarandache Integer Sequences which first appeared in 
Properties of the Numbers, F. Smarandache, University of Craiova Archives, 1975. The first four 
sequences are recurrence generated sequences while the last three are concatenation sequences. 

The Non-Arithmetic Progression: {a; :a; is the smallest integer such that a>a;, and such that 
for ksi there are at most t-1 equal differences a, —@, =@, -@, =...=@,.— A, } 


A strategy for building a t-term non-arithmetic progression is developed and computer implemented 
for 3<t<15 to find the first 100 terms. Results are given in tables and graphs together with some 
observations on the behaviour of these sequences. 

The prime-Product Sequence: {ty : ta = Pr#tl, Pa is the nth prime number}, where p,# 
denotes the product of all prime numbers which are less than or equal to pz. 

The number of primes q among the first 200 terms of the prime-product sequence is given by 6<q<9. 
The six confirmed primes are terms numero 1, 2, 3, 4, 5 and 11. The three terms which are either 
primes or pseudo primes (according to Fermat’s little theorem) are terms numero 75, 171 and 172. 
The latter two are the terms 1019#+1 and 1021#+1. 

: The Square-Product Sequence: {tq! ty = (a!)*+1} 
As in the previous sequence the number of primes in the sequence is of particular interest. Complete 

“prime factorization was carried out for the first 37 terms and the number of prime factors f was 
recorded. Terms 38 and 39 are composite but were not completely factorized. Complete factorization 
was obtained for term no 40. The terms of this sequence are in general much more time consuming to 
factorize than those of the prime-product sequence which accounts for the more limited results. Using 
the same method as for the prime-product sequence the terms t, in the interval 40<n<200 which may 
possible be primes were identified. There are only two of them, term #65: N=(65!)°+1 which is a 182 
digit number and term #76: N=(76!)*+1 which has 223 digits. 

The Prime-Digital Sub-Sequence: The prime-digital sub-sequence is the set 
{M=aotay-10-+a2-10*+...a,10* :M is a prime and all digits ao, a), a2...a, are primes} 

A proof is given for the theorem: The Smarandache prime-digital sub sequence _is infinite, which 
until now has been a conjecture. 

Smarandache Concatenated Sequences: Let G={g1, 82, .-.. Sx, .-.. } be an ordered set of 
positive integers with a given property G. The corresponding concatenated S.G sequence is defined 
through S.G = {a,:a, = g,,4, =@,_,-10°# +9, k> 1}. 

The S.Odd Sequence: Fermat’s little theorem was used to find all primes/pseudo-primes 
among the first 200 terms. There are only five cases which all were confirmed to be primes using the 
elliptic curve prime factorization program, the largest being term 49: 

1357911131517192 1232527293 13335373941434547495 1535557596 1636567697 1737577798 1838587899 1939597 
Term #201 is a 548 digit number. 

The S.Even Sequence: The question how many terms are nth powers of a positive integer was 
investigated. It was found that there is not even a perfect square among the first 200 terms of the 
sequence. Are there terms in this sequence which are 2-p where p is a prime (or pseudo prime)? 
Strangely enough not a single term was found to be of the form 2-p. 

The S.Prime Sequence: How many are primes? Again we apply the method of finding the 
number of primes/pseudo primes among the first 200 terms. Terms #2 and #4 are primes, namely 23 
and 2357. There are only two other cases which are not proved to be composite numbers: term #128 
which is a 355 digit number and term #174 which is a 499 digit number. 


J. The Non-Arithmetic Progression 


This integer sequence was defined in simple terms in the February 1997 issue of Personal Computer 
World. It originates from the collection of Smarandache Notions. We consider an ascending sequence 
of positive integers a), a2, ... a, Such that each element is as small as possible and no t-term arithmetic 
progression is in the sequence. In order to attack the problem of building such sequences we need a 
more operational definition. 


Definition: The t-term non-arithmetic progression is defined as the set : 

{a; :a; is the smallest integer such that a>a,, and such that for k<i there are at most t-1 equal 

differences 2, -a, =a, ~&, =.= G,  - a} 

From this definition we can easily formulate the starting set of a t-term non-arithmetic progression: 
{1, 2,3 ..... t-l, t+1} or {a;: a=i for ist-1 and a=t+1 where t23} 


‘It may seem clumsy to bother to express these simple definitions in stringent terms but it is in fact 
' absolutely necessary in order to formulate a computer algorithm to generate the terms of these 
sequences. 


Question: How does the density of a t-term non arithmetic progression vary with t. i.e. how does the 
fraction a,/k behave for t>3?! 


Strategy for building a t-term non-arithmetic progression: Given the terms a), a2, ... a we will 
examine in turn the following candidates for the term ays: 


ant = a +d, d=1, 2, 3, ... 
Our solution is the smallest d for which none of the sets 

{@), a2, ... a, agtd, a-td-e, a,td-2e, ... a,td-(t-1)-e :e2d} 
contains a t-term arithmetic progression. 
We are certain that a,,; exists because in the worst case we. may have to continue constructing sets 
until the term a;,+d-(t-1)e is less than 1 in which case all possibilities have been tried with no t terms 
in arithmetic progression. The method is illustrated with an example in diagram 1. 
In the computer application of the above method the known terms of a no t-term arithmetic 
progression were stored in an array. The trial terms were in each case added to this array. In the 
example we have for d=1, e=1 the array: 1,2,3,5,6,8,9,10,11,10,9,8. The terms are arranged in 
ascending order: 1,2,3,5,6,8,8,9,9,10,10,11. Three terms 8,9 and 10 are duplicated and 11 therefore 


has to be rejected. For d=3, e=3 we have 1,2,3,5,6,8,9,10,13,10,7,4 or in ascending order: 
1,2,3,4,5,6,7,8,9,10,10,13 this is acceptable but we have to check for all values of e that produce terms 


' This question is slightly different from the one posed in the Personal Computer World where also a wider definition ofa t-term non 
arithmetic progression is used in that it allows a2>a, to be chosen arbitrarily. 
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which may form a 4-term arithmetic progression and as we can see from diagram | this happens for 
d=3, e=4, so 13 has to be rejected. However, for d=5, e=5 no 4-term arithmetic progression is formed 
and e=6 does not produce terms that need to be checked, hence ay = 15. 


Diagram 1.To find the 9” term of the 4-term non-arithmetic progression. 


Routines for ordering an array in ascending order and checking for duplication of terms were 
included in a QBASIC program to implement the above strategy. 


Diagram 2. ax/k for non-arithmetic progressions with !=3, 4, 5, ... 15. Bars are shown for k = multiples of 10. 


Results and observations: Calculations were carried out for 3<t<15 to find the first 100 terms of each 
sequence. The first 65 terms and the 100” term are shown in table 1. In diagram 2 the fractions a,/k 
has been chosen as a measure of the density of these sequences. The looser the terms are packed the 
larger is a,/k. In fact for t>100 the value of a,/k = 1 for the first 100 terms. 


In table 1 there is an interesting leap for t=3 between the 64" and the 65” terms in that ag, = 365 and 
ags = 730. Looking a little closer at such leaps we find that: 
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Table 1. The 45 first terms of the nor-arithmetic progressions for t=3 to 15. 
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Leap starts at Leap finishes at 


5 10 

14 =3-5-1 28 =2-14 
4) =3-14-1 82 =2-41 
122 =3-41-1 244 =2-122 
365 =3-122-1 730 =2-365 


Does this chain of regularity continue indefinitely? 


Sometimes it is easier to look at what is missing than to look at what we have. Here are some 
observations on the only excluded integers when forming the first 100 terms for t=11, 12, 13 and 14. 


For t=11: 11, 22, 33, 44, 55, 66, 77, 88, 99 The n¢h missing integer is 11-n 
For t=12: 12, 23, 34, 45, 56, 67, 78, 89, 100 The nth missing integer is 11-n+1 
For t=13: 13, 26, 39, 52, 65, 78, 91, 104 The nth missing integer is 13-n 
For t=14: 14, 27, 40, 53, 66, 79, 92, 105 The nth missing integer is 13-n+1 


Do these regularities of missing integers continue indefinitely? What about similar observations for 
other values of t? 


IL. The Prime-Product Sequence 


The prime-product sequence originates from Smarandache Notions. It was presented to readers of the 
Personal Computer World’s Numbers Count Column in February 1997. 


Definition: The terms of the prime-product sequence are defined through {t, : t, = Pa#t1, pa is the nth 
prime number}, where p,# denotes the product of all prime numbers which are less than or equal to 


Pn- 


The sequence begins {3, 7, 31, 211, 2311, 30031, ... }. In the initial definition of this sequence t; was 
defined to be equal to 2. However, there seems to be no reason for this exception. 


Question: How many members of this sequence are prime numbers? 


The question is in the same category as questions like ‘How many prime twins are there?, How many 
Carmichael numbers are there?, etc.’ So we may have to contend ourselves by finding how 
frequently we find prime numbers when examining a fairly large number of terms of this sequence. 


From the definition it is clear that the smallest prime number which divides t, is larger than pr. The 
terms of this sequence grow rapidly. The prime number functions prmdiv(n) and nxtprm(n) built into 
the Ubasic programming language were used to construct a prime factorization program for n<10". 
This program was used to factorize the 18 first terms of the sequence. An elliptic curve factorization 
program, ECM.UB, conceived by Y. Kida was adapted to generate and factorize further terms up to 
and including the 49zh term. The result is shown in table 2. All terms analysed were found to be 
square free. A scatter diagram, Diagram 3, illustrates how many prime factors there are in each term . 


The 50th term presented a problem. tso=126173-n, where n has at least two factors. At this point 
prime factorization begins to be too time consuming and after a few more terms the numbers will be 
too large to handle with the above mentioned program. To obtain more information the method of 
factorizing was given up in favor of using Fermat’s theorem to eliminate terms which are definitely 
not prime numbers. We recall Fermat’s little theorem: 
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If p is a prime number and (a, p)=1 then a?' = 1 (mod p). 


a™' = 1 (mod n) is therefore a necessary but not sufficient condition for n to be a prime number. If n 
fills the congruence without being a prime number then n is called a pseudo prime to the base a, 
psp(a). We will proceed to find all terms in the sequence which fill the congruence 


a’"' = 1(modt,) 
for 50 <n < 200. ty is a 513 digit number so we need to reduce the powers of a to the modulus t, 
gradually as we go along. For this purpose we write t,-1 to the base 2: 


trl = > 5(k)-2* , where 5(k) € {0,1} 
k=1 
From this we have 


” k 
qa! = I] qitt)2 


=1 
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Diagram 3. The number of prime factors in the first 49 terms of the prime-product sequence. 


This product expression for a'~' is used in the following Ubasic program to carry out the reduction 


of a"! modulus t,. Terms for which 5(k)=0 are ignored in the expansion were the exponents k are 
contained in the array E%(). The residue modulus t, is stored in F. In the program below the 


reduction is done to base A=7. 


100 dim E%(1000) 

110 M=N-11%=0 

120 T=1J%=0 

130 while (M-T)>=0 

140 ine J%T=2°T 

150 wend 

160 dec J%:M=M-T\2:inc 1%:E%(I%)=J% 
170 if M>0 then goto 120 
180 F=1 

190 for J%=1 to 1% 

200 A=7 

210 for K%=1 to E%(J%) 
240 A=(AA2)QN 

250 next 

260 F=F*A:F=FQN 

270 next 
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Table 2. Prime factorization of prime-product terms 


223092871 =317 - 703763 
6469693231 =331 - 57] - 34231 
200560490131 Prime number 
7420738134811 =181- 60611 - 676421 
304250263527211 =61 - 450451 - 11072701 
130827 61331 670031= 61 - 450451 - 11072701 . 
$3 4889782588491 411 = 953 - 44727 - 13808181181 
325891 58477190044731 = 73 - 139 - 173 - 18564761860301 
19227603501 5421 2639071 =277 - 3467 - 105229 - 19026377261 
117288381359406970983271 =223 - 5259568467082542470777 
7858321 551080267055879091 =547307 29297 - 143581524529603 
S579 408301 266989 609 67 415391 = 1063 - 303049 - 598841 - 2897214489673 
40729 5805992490241 50621323471 = 2521 - 161561 604915704181 4780695! 
32176447 6734067 2907 899084554131 = 22093 - 15031819651 - 9688841 4202798247 
26706451 5689 2758513556240 7992791 = 265739 - 1004988035964897329 I 67431269 
237 $87 41 896345550770650537 601358311 = 131 - 1039 - 2719 - 6422589 1884294373371 806141 
23055579639 4551 84247531021 47331756071 = 2336993 - 13848803 - 7123743602409 1007473549 
232862364358497 36090006331 6880507343071 = 960703 - 242387 46455303809907959 41 27301057 
239848235289 252281 72706521 6386922583962! 1 = 2297 - 9700398839 - 179365737007 - 6001315443334531 
255537 4117594999 4144795978} 534007! 648394471 = 149 - 13203797 - 30501 264491063137 - 42767843651083711 
27973499 681785493651 7827 41 5187 2067809 674997231 = 334507 - 1290433 - 6480464442342997 1 4623377554034701 
3161005454041 76077881 4520629 154346249327 468699 1 = 5122427 - 2025436786007- 3046707595069540247 157055819 
401 447 49393330361 8909-4441 1990260451 3664588524773) = 
1843 - 49999 - 552001 - 57900988201093 -3280805299990739 67231 
52589479052627 7 4077 137179707241 1912700610967 452631 = 
1951 - 22993 - 1172323185947301 41 449323454664 151 43728266617 
720478176302100004856779341 989 20432067 38370254101031 1 = 
881 - 1657 - 32633677 - 160823938621 - 5330079340103 - 1764291759303233 
1001 4646650599 1 90067 509 2331316499 4005735633465320043309 I = 
8782799 5900552888249868 1 487 - 1476476861 45442451 39224580493 
1492182350939 27932005887 573661 5841068547 583863326864530411 = 
87549524399 - 6501816157 3521013453 '> 2621 40076844134219184937113 
22531953499 1831 17732889023622899 2001 3506851 336235654407 1911 = 
23269086799 180847 - 9683213481319911991 636641 541802024271084713 
1353751 669937 174948405357 67087951 74421 2057570647 889977 42242987 1 = 
1381 - 1867 - 8311930927 - 38893867968570583 - 42440201875440880489 1 13304753 
576615221997 5951 659023630035336 1 343065653840 I 560606631 985606881 1 = 
1361 - 214114727210560829 - 32267019267402210517 - 613228865630544238382107 
952947 4207 359839270569 4621 5901 1344291954191 3060621 3075415763491 271 = 
205590139 - 53252429177 - 70645763395667 63 - 124501547099289 40906 1 979 46067 239 
16658990378732521 938085 1 69535089 6256250980509 59 487 48620469 61 683989711 = 
52614127 - 26605801 5409361 15803523331 939275661 5852809877 22606890621 81793 
2981 959277793121 4269 1724534678104298689 2551 12174826003064061 41434158091 = 
$01 - 1653781 - 8564177 - 358995947 - 1525310189119 - 640532866409 661 895480902986 1252251 
539734629 28055497827 2021 4077 673687 80627 5517 530364350655459 51 1599582614291 = 
107453 - 5634838141 - 891 4157280964101 12334489 1396571 2571 6363297 4628403 1 7 4028667 
10308931 41925860008499 56088383567 437099862384829959097 SI 927 667 1 SS20279329391 = 
32999 - 175603474759 - 771 48541513247 - 2305961 4664373239 5959853041 586242331 6227152033 
19896237639 1 690981 640415251 5452851 536027344027 2182105821 220397 6095413910572271 = 
21639496447 - 79791 25905957339 49501 8877 - 1152307771 6259797580440201 62101777 453615909 
391955881491 631 233831 51804554421 175259738677 3361 987 48447 8041 8329079 6540382737191 = 
$21831 - 50257723 - 1601684368321 - 39081 170243262541027 - 23875913958369977 1585726531 60969521 
77999 22041 683451 553249199 3063298 | 387 668799 67899035509 4509303247 486851 1536164700811 = 
467 - 10723 - 57622771 » 5876645549 - 9458145520867 - 486325954430626096097 1 9222040521 49 47865503847 
164578355079521038773558101 1 435590727981 167322669 64924941 4629852197 255934130751870911 = 
1051 - 2179 - 14333 - 43283699 - 75311908487 - 29281 2710684839 - 4609659657 286645929343033404487 2907 384889 
3670097318273319 1 6465034565550! 3467323398003 1 295533178261 946245703998807 331 1157667212931 = 
13867889 4681 59 - 2646471 42357 1660867 5791 59849289 6703564888 1 000360533429306 1 9468037 57 2880509 
83311209 124804345037 56284637988 1038241 13467104086031 46546179777 48077 292641 632790457335 111 = 
3187 . 31223 - 1737142793 - 1146303934031 54601 - 973104505470446969309 113 - 43206785807 56718923287 5099500379 
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This program revealed that there are at most three terms t, of the sequence in the interval 50<n<200 
which could be prime numbers. These are: 


Term #75. N=379#+1_N is a 154 digit number. 
N=171962010545840643348334056831 754301 95845756358957 42560438771 10505832! 6552385626 1308397965147 
955578800999 45578220245652269329062952082627568222756636941 11 
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Term #171, N=1019#+1.N is go 425 digit number. 

N=2040406899301 6374194542464] 727746076956597971 174231 219132271 310323390261 69 1 759299022444537574 
1046872884292986227 1 60555781 8821 68549067666 1 9853898399 58622802465986881 3761 394041383761530961031 
408346655636467 401 6027975521 231 7501 3568630036386 I 239066 1 66840623542231 17837423905 105265872570265 
0030269683479324852673430580 1 6341 65948702506367 176701 23329806461 66635537 1 697542904875 15755971504 
17381063934255689 | 24486029 4929089 66644747931 


Term #172. N=1021 #+1.N is a 428 digit number. 

N=208325544418697 1805262785592040287 4457 268652856889007 47340490078401 81 457187286244301 915872863 
160885721 486313893793092847 4301 694088598087 1 887083026597 75388131 77 7260588503833 1 6252820523111 213 
0679219354048332 1 7036456300717761 68885357 1 267 15023250865563442766366 1 80331 200980711 2476455894240 
568090534683239067 45795726223468483433625259000887 41 1959 1973239736 1348834503191 30587753586846905 
7614606627 6875058596 1002361 12260054944287636531 


The last two primes or pseudo primes are remarkable in that they are generated by the prime twins 
1019 and 1021. 


Summary of results: The number of primes q among the first 200 terms of the prime-product 
sequence is given by 6<q<9. The six confirmed primes are terms numero 1, 2, 3, 4, 5 and 11. The 
three terms which are either primes or pseudo primes are terms numero 75, 171 and 172. The latter 
two are the terms 1019#+1 and 1021#+1. 


YL The Square-Product Sequence 


Definition: The terms of the square-product sequence are defined through {t,: ta = (n!)°+1} 


This sequence has a structure which is similar to the prime-product sequence. The analysis is 
therefore carried out almost identically to the one done for the prime-product sequence. We merely 
have to state the results and compare them. 


The sequence begins {2, 5, 37, 577, 14401, 518401, ... } 


As for the prime-product sequence the question of how many are prime numbers has been raised and 
we may never know. There are similarities between these two sequences. There are quite a few primes 
among the first terms. After that they become more and more rare. Complete factorization of the 37 
first terms of the square-product sequence was obtained and has been used in diagram 4 which should 
be compared with the corresponding diagram 3 for the prime-product sequence. 
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Diagram 4. The number of prime factors in the first 40 terms of the square-product sequence. 
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Diagram 4 is based on table 3 which shows the prime factorization of the 40 first terms in the square- 
product sequence. The number of factors of each term is denoted f. The factorization is not complete 
for terms numero 38 and 39. A +-sign in the column for f indicates that the last factor is not a prime. 
The terms of this sequence are in general much more time consuming to factorize than those of the 
prime-product sequence which accounts for the more limited results in this section. Using the same 
method as for the prime-product sequence the terms t, in the interval 40<n<200 which may possible 
be primes were identified. There are only two of them: 


Term #65. N={(65!)2+1. N is a 182 digit number. 
680237 40289078328950450781 97 262220379290257 695327 | 358034279380 1 04027 100652464382649659623724446578151 4128589 


9657 153438534056379 295 | 822384455 | 8074780057 600000000000000000000000000000 1 


Term #76 N=(76l)2+1.N is a 223 digit number. 
35550902700 107478542025131357707726481 9432566692554 1 64797 7005 2502800500841 772266884421 39 166589065 1 6439209 129 


30369944999 45253 100626495077 57826978507 1 9865801 162529840993 1764786386381 1506 1 7600000000000000000000000000000 
0000001 


Table 3. Prime factorization of square-product terms. 


Ha oe Re[ni}#1 and fe factor 
1 [2 
5 
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1625702401= 17- 95629553 
131681894401 
13168189440001 
1593350922240001 
229442532802560001=101- 2271708245569901 
38775788043432640001 
7600054456551997 440001= 29- 109- 2404319663572286441 
17100122527241 99424000001 =13441 69-1 2721 70577304043929 
4377631 3669739 5052544000001 =1 49- 2938007628841 577533852349 
126513546505547 17018521 6000001= 9049-13980942259 4261 4324071 3449 
40990389057797 2831 4000998400000 1 = 37-1 1078483531 83710355135404972973 
1479753045347 481 9213543604224000001=7103-41 -2083 I $87 1 58104092560535861 261 
5919012181389927 685417441 689600000001=41-10657- 8681601 7- 348046955409-4483247 49841 
261028437 1992958 10926909 17851 1360000000!= 63-157-272557 6247 25170524096177 486176631513 
1263377 63604459 172488624042399 4982400000001 = 337- 8017- 51 4049836440277 481- 90967 4823323537849 
$683267 69467589022 464821 184293345589 600000001 = 509-1 544837 4629- 8499 40026045327 47687 401741723441 
3849562192133312769397370021 529 67 1 1720960000000 1 
240597 637008332048087323562634560444825600000000000 I = 
941- 813769831 908479758733- 313425331349331 290243399417 
162644002617632464507038883409 62860702 I 0S600000000000 I= 
53- S3- 418633- 6017159668589- 229858897 1 287 609 6222556462301797 
118567477908254065625631 3460056192545 1834982400000000000 | = 
113 42461- 745837. 2460281- 7566641- 15238649- 116793504008451 126962009 
9295690268007 | 18823449 497526840549 554238626201 00000000000 I = 
212259034557 6634509- 43794085292997 939303952241 47 49827 53464389 
781767551 5393986930521027 4200729021751 14684635545600000000000 1 = 
1713 7078604732075883.49837- 455289320701 41406371 6469396531 758248773 
70359079 53854588237 46892457 806561 1957460321 6171991 040000000000000 I = 
$1- 1733 15661- 359525849- 100636381 1265686901 10069- 1174592249518207759S37897 
$761 507 55326425929 6207 634615621 05309 1 2566907 41 283389 440000000000000 I = 
353- 422041- 13400767181 - 338676081 809 484090853058207938321 9 1570324667821 677 
$923783734542601 51931 661 989439595834544685 13190741907 86560000000000000 1 = 
10591621 681- 6415450838021. 522303293914660001204969- 1950882388585355532025429 
75400004869 1 6893054535799064997 19865997 1 62108647 17937 665$3840000000000000 1 = 
37- 3121- 4421- 4073332882845935253- 362581351 23244480427 3874507 621085782231 48052301 
871 624056287 592837 10433837 191 36751 650927 193975961393594 1 47799040000000000000 I = 
193- 13217- 866100731693- 394521434431 45645231 47689 4644096901 29 1 1974106242868 1 6576197 
1067739 468952301 22545281 450559 425330223858 1 2620552707 152831 0538240000000000000001 = 
317 - 373 - 9030196538846808-4889 78 28545563235536086347 48201 1761 632198907771 618981 5715341 
1383790351762182388 1868475992501 52279701 20131562363084700690457 55904000000000000000 I = 
73- 57986941373 - 3269017431 698277804505207 628624951781 7882480631366! 4675401 16377341 49869 
189440899 1 562.427 689 42779 43633734584709 109 44601 0887 $0629 552452363983257 5000000000000000 I = 
127 404354297881 - 49105571 194338128021910109 - 302797201 1 452403842829 24307698 | 4272052327 543069 
273552658381 $1 455835337350607 1 1 27 4031 99540400397 21 5590907 37 41 21 359 182397 44000000000000000 I = 
233-757- 15509190807 49 1428121 7009 4885906800637 254241 6722731 7903236388341 9 1845062531 6785821 6021 
41 60735933984357 43255481 102734 1847802665009 49004 | 649 137701 1 603858731 642650624000000000000000 I = 
61-10045457577 41-679001 282693270665470845882 1 250934897329 61 $34897 2937 285510229243 1 $7 27303361 7801 
5657177 49437 497 1892087 5976437 469 5648426401 51 8406663862032 185661 7397062824099 8400000000000000000 I = 
89- 701- 187100101949. 570306192879849 15195631 567 31 422223621 393496539544379 4926477 2817 4834902025511 
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IV. The Smarandache Prime-Digital Sub-Sequence 


Definition: The prime-digital sub-sequence is the set {M=aota;-10+a2- 107+...a,10* :M is a prime and 
all digits ao, a), a2...d, are primes} 


The first terms of this sequence are {2, 3, 5, 7, 23, 37, 53, 73, ... }. Sylvester Smith [1] conjectured 
that this sequence is infinite. In this paper we will prove that this sequence is in fact infinite. Let’s 
first calculate some more terms of the sequence and at the same time find how many terms there is in 


the sequence in a given interval, say between 10* and 1 


om! 


.The program below is written in Ubasic. 


One version of the program has been used to produce table 4 showing the first 100 terms of the 
sequence. The output of the actual version has been used to produce the calculated part of table 5 


which we are going to compare with the theoretically estimated part in the same table. 
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Ubasic program 
10 point 2 
20 dim A%(4),B%(4) 
30 for 1%=1 to éread AR%(I%):next 
40 data 1,4,6,8,9,0 ‘Digits not allowed stored in A%() 
50 for I%=1 to 4read B%(I%):next 
60 data 23.5.7 ‘Digits allowed stored in B%() 
70 for K%=1 to7 ‘Calc. for 7 separate intervals 
80 M%=0:N=0 
90 for E%=1 to 4 ‘Only 2,3,5 and 7 allowed as first digit 
100 P=B%(E%)*10AK%:PO=P:S=(B%[E%) +1 )*]OAK%:gosub 150 
110 next 
120 print K%,.M%,N,M%/N 
130 next 
140 end 
150 while P<S 
160 P=nxtporm(P):P$=sir(P) ‘Select prime and convert to string 
170 incN ‘Count number of primes 
180 L%=len{P$):C%=0 ‘C% will be set to 1 if P nof member 
190 forlZ=2tol% 
200 for J%=1to 6 ‘This loop examines each digit of P 
210 if val{mid(P$,1%,1))=A%(J%) then C%=1 
220 nextnext 
230 if C%=O then inc M% ‘If criteria filed count member (m%) 
240 wend 
250 retum 
Table 4. The first 100 terms in the prime-digital sub sequence. 
a 
2 3 5 7 23 37 53 
233 257 277 337 353 373 $23 
733 757 773 2237 2273 2333 2357 2377 


Table 5. Comparison of results. 


4 38 389 1325 

0.6021 1.1761 1.5798 2.1072 2.5899 3.1222 
13 472 3771 30848 261682 
0.30769 0.08051 0.03394 0.01261 0.00506 


34 364 
1.5278 2.5615 3.6430 
421 28464 2146190 
0.08000 0.01280 0.00205 


Theorem: 
The Smarandache prime-digital sub sequence is infinite. 


Proof: 
We recall the prime counting function m(x). The number of primes p<x is denoted (x). For 


— x 
sufficiently large values of x the order of magnitude of 7(x) is given by 2(x) = ieee: Let a 


and b be digits such that a>b+0 and n(a,b,k) be the approximate number of primes in the interval 
(b-10*,a-10*). Applying the prime number counting theorem we then have: 


Ww _t a () 
I 7 logd 
x logl0+—>= logl0-+—— 


n(a,b,k) = 


Potential candidates for members of the prime-digital sub sequence will have first digits 2,3,5 or 7, ie. 
for a given k they will be found in the intervals (2-10*,4-10*), (5-10*,6-10*) and (7-10*8-10"). The 
approximate number of primes n(k) in the interval (10*,10") which might be members of the 
sequence is therefore: 


n(k)=n(4,2,k)+n(6,5,k)+n(8,7,k) (2) 


The theoretical estimates of n in table 5 are calculated using (2) ignoring the fact that results may not 
be all that good for small values of k. 


We will now find an estimate for the number of candidates m(k) which qualify as members of the 
sequence. The final digit of a prime number >5 can only be 1,3,7 or 9. Assuming that these will occur 
with equal probability only half of the candidates will qualify. The first digit is already fixed by our 
selection of intervals. For the remaining k-1 digits we have ten possibilities, namely 0,1,2,3,4,5,6,7,8 
and 9 of which only 2,3,5 and 7 are good. The probability that all k-1 digits are good is therefore 
(4/10)". The probability q that a candidate qualifies as a member of the sequence is 
re eee 

q=5°Gp (3) 
The estimated number of members of the sequence in the interval (10*,10"*') is therefore given by 
m(k)=q-n(k). The estimated values are given in table 5. A comparison between the computer count 


and the theoretically estimated values shows a very close fit as can be seen from diagram 5 where 
logio m is plotted against k. 


The prime-digital sub sequence 


Diagram 5. logiom as a function of k. The upper curve corresponds to the computer count. 


loga 
k 


For large values of k we can ignore the terms 


For large k we therefore have 


(a —6)10* 
mA DE) klogl0 


and (2) becomes 


4-10* 


ha 
me) klog10 


Combining this with (3) we get 


ce Na 


2 Saar SPST 


d == in comparison with log 10 in (1). 
(1’) 
(2’) 
(4) 


From which we see (apply for instance |’Hospital’s rule) that m(k)—><0 as k->eo. A fortiori the prime- 


digital sub sequence is infinite. 


V. Smarandache Concatenated Sequences 


Smarandache formulated a series of very artificially conceived sequences through concatenation. The 
sequences studied below are special cases of the Smarandache Concatenated S-sequence. 


Definition: Let G={g,, g2, .... 2x .... } be an ordered set of positive integers with a given property G. 
The corresponding concatenated S.G sequence is defined through 
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S.G = {a,:4, = g,,a€, =a,_,-10°# +9, k>D 


In table 6 the first 20 terms are listed for three cases, which we will deal with in some detail below. 


Table 6. The first 20 terms of three concatenated sequences 


The S.odd sequence 


1 2 2 
24 


246 

2463 

246810 

24681012 
135791313 2468101214 2357111317 
13579111315 2468101 21416 Z3S7VVISIT7I9 
135791 1131517 2468101 2141618 23571113171923 
135791 113151719 246810121 4161820 2357111317192329 


235711131719232931 
23571133171923293137 
2357111317192329313741 
235711131719232931374143 
2357111317192329313741 4347 
2357111337392329313741 434753 

23571 11317192329313741 43475359 
2357111317192329313741 4347535961 
2357111317192329313741 434753596167 
2357111317192329313741 4347535961 6771 
2357111317192329313741 4347535961677173 


246810121 416182022 

246810125 41618202224 

246810121 4161820222426 

246810321 43 6182022242628 

246810121 41 618202224262830 

2468101 21 4161820222426283032 
246810121 416182022242628303234 
2468103 21 41 618202224262830323436 
2468101 21 4161820222426283032343633 
246810121 41 6182022242628303234363840 
246810121 4161 8202224262830323436384042 


13579111315171921 

335791 1131517192123 
135791133351719212325 

135791 13315171921 232527 

135793 1131517192123252729 

135791 113151719212325272931 

135791 11315171921 232527293133 
135791 11315171921 2325272931333S 
13579111315171921 232527 2931333537 
135791 11315171923 2325272931 33353739 
135791 1131517192123252729313335373941 


Case I. The S.odd sequence is generated by choosing G={1,3,5,7,9,11,.....}. Smarandache asks how 
many terms in this sequence are primes and as is often the case we have no answer. But for this and 
the other concatenated sequences we can take a look at a fairly large number of terms and see how 
frequently we find primes or potential primes. As in the case of prime-product sequence we will resort 
to Fermat’s little theorem to find all primes/pseudo-primes among the first 200 terms. If they are not 
too big wee can then proceed to test if they are primes. For the S.odd sequence there are only five 
cases which all were confirmed to be primes using the elliptic curve prime factorization program. In 
table 7 # is the term number, L is the number of digits of N and N is a prime number member of the 
S.odd sequence.: 


Table 7. Prime numbers in the S.odd sequence 


# L N 

2 2 13 

10 1S T3S791HISISI719 

16 27 135791113151719212325272931 

34 6&3 13579111315171921232527293 1333537394 14345474951 5355575961 636567 

49 93 1357911 13151719212325272931 333537394 14345474951 5355575961 636567697 1737577798 1838587899 1939597 


Term #201 is a 548 digit number. 


Case 2. The S.even sequence is generated by choosing G={2,4,6,8,10, ...... }. The question here is : 
How many terms are nth powers of a positive integer? 


A term which is a nth power must be of the form 2"-a where a is an odd nth power. The first step is 
therefore to find the highest power of 2 which divides a given member of the sequence, i.e. to 
determine n and at the same time we will find a. We then have to test if a is a nth power. The Ubasic 
program below has been implemented for the first 200 terms of the sequence. No nth powers were 
fond. 


Ubasic program: (only the essential part of the program is listed) 


60 N=2 

70 for U%=4 to 400 step 2 

80 D%=int(log(U%)/log{10))+1 ‘Determine length of U% 
90 N=N*1OAD%+U% ‘Concatenate U% 

100 A=N:E%=0 

110 repeat 

120 Al=A:A=A\2:inc E% ‘Determine E% {=n) 

130 until res<>0 


? 


132 dec E%:A=Al ‘Determine A (=a) 

140 B=round{AA(1/E%)) 

150. if BAE%=A then print E%,N ‘Check if ais a n/Apower 
160 next 

170 end 


So there is not even a perfect square among the first 200 terms of the S.even sequence. Are there 
terms in this sequence which are 2-p where p is a prime (or pseudo prime). With a small change in the 
program used for the S.odd sequence we can easily find out. Strangely enough not a single term was 
found to be of the form 2-p. 


Case 3. The S.prime sequence is generated by {2,3,5,7,11, ...}. Again we ask: - How many are 
primes? - and again we apply the method of finding the number of primes/pseudo primes among the 
first 200 terms. 


There are only 4 cases to consider: Terms #2 and #4 are primes, namely 23 and 2357. The other two 
cases are: term #128 which is a 355 digit number and term #174 which is a 499 digit number. 


#128 

2357111317192329313741 434753596 16771 7379838997 101 1031071091 13127131 1371391 4915115716316717317918 
119119319719921 1 223227229233239241 25 125726326927 1 277281 28329330731 131331 733133734734935335936737 
337938338939 7 401 409419421 431 433439443449 457 46 1 463467479 487 49 1 499503509521 523541 547557563569571 57 
7587593599601 60761361 7619631641 643647653659661 67367768369 I 701709719 


#174 

2357111317192329313741 4347535961577 1 7379838997101 1031071091 13127131 137139149151 15716316717317918 
1191 19319719921 1223227229233239241 25 1 25726326927 1277281 28329330731 131331 733133734734935335936737 
3379383389397 401 409419421 431 433439443449457 46 | 463467479487 49 | 499503509521 523541 54755756356957157 
7587593599601 60761361 7619631641 6436476535966 1 67367768369 I 701 7097 19727733739 743751757761 76977378 
779780981 182182382782983985385 7859863877881 8838879079 } 1919929937941 947953967971 97798399 | 99710091 
01310191021 10311033 


Are these two numbers prime numbers? 


